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ABSTRACT 

An extensive computational study of supersonic quasi- 
axisymmetric vortex breakdown in a configured circu- 
lar duct is presented. The unsteady, compressible, full 
Navier-Stokes (NS) equations are used for the present 
computational study. The NS equations are solved for 
quasi-axisymmetric flows using an implicit, upwind, flux- 
difference splitting, finite-volume scheme. The quasi- 
axisymmetric solutions are time accurate and are obtained 
by forcing the components of the flowfield vector to be 
equal on two axial planes, which are in close proximity of 
each other. The computational study addresses the effect 
of the Reynolds number, for laminar flows, on the evo- 
lution and persistence of vortex breakdown. The effect 
of boundary conditions at the duct exit on vortex break- 
down is also studied. Finally, the effect of swirl ratio at 
the duct inlet is investigated. 

Introduction 

The majority of the experimental studies of vortex 
breakdown phenomenon has been focused on incompress- 
ible flows in pipes 1 " 5 . Two main types of vortex break- 
down have been documented experimentally: the quasi- 
axisymmetric bubble type and the asymmetric spiral type. 
Other types of vortex breakdown were also generated in 
pipes. 

The major effort of numerical study of vortex break- 
down flows, has also been focused on incompressible, 
quasi-axisymmetric isolated vortices. Grabowski and 
Berger 6 were the first to use the incompressible, quasi- 
axisymmetric NS equations to study isolated vortex flow 
in an. unbounded region. Hafez, et al 7 solved the incom- 
pressible, steady, quasi-axisymmetric Euler and NS equa- 
tions using the stream function-vorticity formulation for 
isolated vortex flows. They predicted vortex breakdown 
flows similar to those of Garbowski and Berger. Salas 
and Kuruvila* solved the unsteady, quasi-axisymmetric 
NS equations in a straight circular pipe. They ob- 
tained steady, multiple bubble-type vortex breakdown for 
a Reynolds number range of 100-1,800. Menne 9 has also 
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used the stream function-vorticity formulation for study- 
ing unsteady, incompressible quasi-axisymmetric isolated 
vortex flows. Wu and Hwang 10 used the stream function- 
vorticity formulation to study quasi-axisymmetric vortex 
breakdown in a pipe. Their study focused on the effects 
of inflow, wall boundary conditions and Reynolds num- 
ber on breakdown structure. They showed that the evo- 
lution of breakdown can be steady, periodic or unsteady 
depending on the inflow velocity profiles and Reynolds 
number. Menne and Liu 11 integrated the laminar, in- 
compressible, NS equations for the breakdown of a vor- 
tex in a slightly diverging pipe. They showed break- 
down flow cases which are based on the purely quasi- 
axisymmetric and non-ax isymmetric analyses. The re- 
sults were in good agreement with the experimental re- 
sults of Leibovich 4 . Spall, et. al 12 used the vorticity- 
velocity formulation of the incompressible NS equations 
to study the three-dimensional vortex breakdown. Breuer 
and Hanel 13 solved the unsteady incompressible NS equa- 
tions using a dual-time stepping, upwind scheme to study 
the temporal evolution of the three-dimensional vortex 
breakdown. In Refs. 12 and 13, both types of breakdown; 
the bubble type and the spiral type, were predicted. Re- 
views of the physical and computational aspects of the in- 
compressible vortex breakdown were presented by Krause 
in Refs. 14 and 15. One of the most important aspect 
of vortex breakdown which Krause discusses in Ref. 15 
is the effect of side boundary conditions on the up- or 
downstream motion of the breakdown point Also, he 
presents alternative outflow boundary conditions. 

Transonic and supersonic flows around highly swept 
wings and slender wing-body configurations at moderate 
to high angles of attack are characterized by interacting 
vortex cores and shock waves. Other applications which 
encounter vortex-shock interaction include a supersonic 
inlet ingesting a vortex and injection of a swirling fuel jet 
into a supersonic combustor to enhance mixing between 
the air stream and the fuel 16 ' 11 . In all these applications 
vortex breakdown due to the interaction of a longitudi- 
nal vortex and a shock wave may occur or is intended 
to occur. For such problems, computational schemes are 
needed to study, predict and control voitex-shock inter- 
action including vortex breakdown. Unfortunately, the 
literature lacks this type of analysis with the exception 
of the preliminary work of Liu, Krause and Menne 19 , 
Copening and Anderson 20 , Delery, et al 16 , Kandil and 
Kandil 21 and Meadows, Kumer and Hussaini 22 . 



The first lime-accurate NS solution for a supersonic 
vortex breakdown was developed by the present au- 
thors in Ref. 23. They considered a supersonic quasi- 
axisymmetric vortex flow in a configured circular duct 
The time-accurate solution of the unsteady, compress- 
ible NS equations was obtained using an implicit, up- 
wind, flux -difference splitting finite-volume scheme. A 
shock wave has been generated near the duct inlet and 
unsteady vortex breakdown has been predicted behind 
the shock. The predicted flow was characterized by the 
evolution, convection and shedding of vortex breakdown 
bubbles. The Euler equations were also used to solve the 
same problem. The Euler solution showed larger size and 
number of vortex-breakdown bubbles in comparison with 
those of the NS solutions. The time-accurate solution was 
carried out for 3,200 time steps which are equivalent to 
a dimensionless time of 16. Only one value of Reynolds 
number of 10,000 was considered in Ref. 23. 

In the present paper, we expand our study of this 
flow using time-accurate computations of the NS equa- 
tions with a fine grid in the shock-vortex interaction re- 
gion and for long computational times. Several issues are 
addressed in the present study. First, we show the effect 
of Reynolds number on the temporal evolution and per- 
sistence of vortex-breakdown bubbles behind the shock. 
In this stage of computations, the conditions at the down- 
stream exit are obtained by extrapolating the components 
of the flowfield vector from the interior cell centers. Al- 
though the flow is supersonic over a large portion of the 
duct exit, subsonic flow exists over a small portion of the 
exit around the duct centerline. Therefore, selected flow 
cases have been recomputed using a Riemann-invarient- 
type boundary conditions at subsonic points of the duct 
exit Finally, the effect of swirl ratio at the duct inlet has 
been investigated. 


Formulation 


The conservative form of the dimensionless, un- 
steady, compressible, full NS equations in terms of time- 
independent, body-conformed coordinates and 
is given by 


dQ dE m 

dt dt™ 
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The first element of the three momentum elements of Eq. 
(5) is given by 
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The second and third elements of the momentum elements 
are obtained by replacing the subscript 1, everywhere in 
Eq. (7), with 2 and 3, respectively. The last element of 
Eq. (5) is given by 
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The reference parameters for the dimensionless form 
of the equations are L, L/a^, poc and for the 
length, velocity, time, density and molecular viscosity, 
respectively. The Reynolds number is defined as Re = 
PocVocL/jioc, where L is the initial radius of the vor- 
tex or the duct inlet radius. The pressure, p, is related 
to the total energy per unit mass and density by the gas 
equation 


P = (T-l)p[e-I(uJ + u| + u^)j (9) 

The viscosity is calculated from the Sutherland law 

/1 =tJ/j (tt§) , c= 04317 (10) 


and the Prandtl number P f = 0.72. In Eqs. (l)-{8), the 
indicial notation is used for convenience. 






Computational Scheme 

The Computational scheme used to solve the un- 
steady, compressible full NS equations is an implicit, 
upwind, flux -difference splitting, finite-volume scheme. 
It employs the flux-difference splitting scheme of Roe 
which is based on the solution of the approximate one- 
dimensional Riemann problem in each of the three di- 
rections. In the Roe scheme, the inviscid flux differ- 
ence at the interface of a computational cell is split into 
left and right flux differences. The splitting is accom- 
plished according to the signs of the eigenvalues of the 
Roe averaged- Jacobian matrix of the inviscid flux at the 
cell interface. The smooth limiter is used to eliminate os- 
cillations in the shock region. The viscous and heat-flux 
terms are linearized and the cross-derivative terms of the 
viscous Jacobians are dropped in the implicit operator. 
These terms are differenced using secord-order spatially 
accurate central differencing. The resulting difference 
equation is approximately factored and is solved in three 
sweeps in the and £ 3 directions. The scheme is 
used for third-order spatial accuracy and first-order tem- 
poral accuracy. The scheme is coded in the computer 
program which is called "FTNS3D”. 

Computational Study and Discussion 

Figure 1 shows a configured circular duct which con- 
sists of a short, straight cylindrical part at the inlet which 
is followed by a short divergent cylindrical part until the 
axial length of 0.74. The divergence angle is 6°. The duct 
radius is then kept constant and a convergent-divergent 
nozzle with a throat radius of 0.95 is attached. The duct 
exit radius is 0.98 and its total length is 2.9. The diver- 
gent part of the duct ensures the stability of the formed 
shock in the inlet region. The configuration of the duct 
is intended to ensure that the supersonic inflow will be- 
come supersonic at the exit As the computations will 
show, a small portion of the duct exit near its center- 
line becomes subsonic at certain times for the specified 
inflow conditions. 

This configured duct has also been used by Delery, et 
al 16 for their Euler equations computations of supersonic 
vortex breakdown in an attempt to computationally model 
an experimental set up. It should be pointed here that 
the Euler equations, used by Delery, et al, assume 
isenthalpic flow in order to drop the energy equation. 
This is a serious approximation since the upstream flow 
is rotational. Moreover, as our present calculations show, 
the flow is actually unsteady and hence, the isenthalpic 
assumption is not valid. 

The NS solver is used to compute all the following 
flow cases by using a grid of 221x51 on two axial planes, 
where 221 points are in the axial direction and 51 points 
are in the radial direction. In the inlet region up to 
the 0.74 axial station, 100 grid points are used and the 
remaining 121 points are used in the remaining part of 
the duct. The grid is also clustered at the centerline (CL) 


and the wall. The minimum radial grid size at the CL is 
0.002. The two axial planes are spaced circumferentially 
at a certain angle so that the aspect ratio of the minimum 
grid size will be less than 2. The upstream Mach number 
is kept at 1.75 and the Reynolds number is varied from 
2,000 up to 100,000. The initial profile for the tangential 
velocity is given by 


w 




1 - exp 
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where Uoc = 1.74, r m = 0.2 and k* = 0.1. The maximum 
swirl ratio /?, is at r = 0.224. The radial velocity, v, 
at* the initial station is set equal to zero and the radial 
momentum equation is integrated to obtain the initial 
pressure profile. Finally, the density p is obtained from 
the definition of the speed of sound for the inlet flow. 
With these compatible set of profiles, the computations 
are carried out accurately in time with At = 0.0025. 
The wall boundary conditions follow the typical Navier- 
Stokes solid-boundary conditions. These computations 
have been carried out on the CRAY YMP of the NASA 
Langley Research Center. The CPU time is 40 ^s/grid 
point/iteration for the NS calculation. 

Next, we present the results of the computational 
study which covers the effects of Reynolds number, the 
exit boundary conditions and the inlet swirl ratio. 


Effect of Reynolds Number 

For these flow cases, the Reynolds number values are 
2,000; 4,000; 10,000; 20,000 and 100,000. The swirl 
ratio, /?, is kept fixed at 0.32. The exit boundary condi- 
tions are obtained by extrapolation from the interior cell 
centers. 

Re = 2,000 

Figure 2 shows the streamlines and Mach contours for 
this flow case at t = 11, which is equivalent to 4,400 
time steps. No vortex breakdown develops and the Mach 
contours show a steady shock at the duct inlet The shock 
is a normal shock over most of the duct inlet The flow 
at the duct exit is supersonic. 

Re = 4,000 

Figure 3 shows snapshots of the streamlines and Mach 
contours for the flow case of Re = 4,000. For this 
value of Reynolds number a single breakdown bubble is 
seen at t = 5 and it is convected downstreams as time 
passes. This breakdown bubble is formed during the 
downstream motion of the inlet shock, which reaches its 
maximum downstream displacement at t = 5. Later on, 
the shock moves upstream, as it is seen at t = 8, while 
the breakdown bubble is convected in the downstream 
direction. Thereafter, the shock stays stationary at the 
inlet This swirling flow case shows a transient single 
breakdown flow. It should be noticed that at t = 5 a 
small portion erf the duct exit at its centerline becomes 
subsonic. At t = 8, it expands radially to about 25% of 
the duct exit radius. 

Re = 10,000 

Figures 4 and 5 show snapshots of the streamlines and 





Mach contours for the flow case of Re » 10,000. At 
t = 3, a single breakdown bubble is formed behind the 
downstream moving shock. In this range of t = 3-5, the 
bubble grows in all directions while the shock moves 
downstreams. In the time range of t = 6-8, the break- 
down bubble splits into two bubbles which are convected 
as they diffuse in the downstream direction. The inlet 
shock moves upstreams during this time range. During 
the time range of t = 9-12, the inlet shock slightly oscil- 
lates at the duct inlet while a new breakdown bubble is 
formed behind the shock. As the bubble size and strength 
increase, the inlet shock moves very slowly downstreams. 
More breakdown bubbles (three bubbles) are formed as 
seen at t = 17. The breakdown bubbles are then con- 
vected downstreams (t = 19) while the inlet shock again 
moves upstreams. As the most downstream bubble is 
shed through the totally supersonic exit, the inlet shock 
slowly moves downstreams. In the time range t = 22- 
33, the shock keeps its slow downstream motion until it 
becomes stationary. During this time, downstream shed- 
ding of breakdown bubbles continues and the formation 
of new breakdown bubbles slows down until it stops when 
the inlet shock becomes stationary. Thereafter, no break- 
down bubbles are formed and the flow becomes steady. 
This swirling flow case shows a transient multi-bubble 
breakdown flow. 

Re * 20,000 

Figures 6 and 7 show snapshots of the streamlines and 
Mach contours for the flow case of Re = 20,000. The 
mechanism of evolution, convection and shedding of the 
vortex-breakdown bubbles while the inlet shock is mov- 
ing downstreams, then upstreams and finally downstreams 
to become stationary is very similar to the previous case 
of Re = 10,000. However there are some few differences. 
First, the size, number and strength of breakdown bubbles 
are larger than those of the case of Re = 10,000. Second, 
the displacements of the inlet shock are larger than those 
of the case of Re = 10,000. Third, the transient time of the 
multi-bubble breakdown is longer than that of the case of 
Re = 10,000. The reader can easily compare the snapshots 
of the streamlines and Mach contours of the two cases at 
exactly the same time instants. Again this swirling flow 
case shows a transient multi-bubble breakdown flow. 

Re = 100,000 

Figures 8 and 9 show snapshots of the streamlines and 
Mach contours for the flow case of Re » 100,000. It 
should be noted that the radial extension of the shown 
streamlines snapshots is r = 0.6, which is larger than 
those of Re = 2,000 - 20,000. Again the mechanism of 
evolution, convection and shedding of the vortex break- 
down bubbles up to t = 30 is very similar to the previous 
flow cases of Re = 10,000 and 20,000. It is noticed that 
the size, number and strength of breakdown bubbles are 
larger than those of the previous cases. Moreover, it is 
noticed here that short periodic evolution, merging, con- 
vection and shedding cycles of the breakdown bubbles 
occur, e.g.; the time periods of 16-21, 22-27 and 28-32. 


At t = 33 and beyond, a new mode of evolution, convec- 
tion and shedding of the breakdown bubbles occurs. It 
should be noticed that the inlet shock keeps on moving 
slowly in the downstream direction and another shock, 
which is downstream of the inlet shock, does not extend 
to the duct axis and is first seen at t = 30, also keeps on 
moving in the downstream direction. During the motion 
of these shocks a breakdown bubble grows behind the 
inlet shock and sheds a breakdown bubble in the down- 
stream direction. At t = 61, the downstream bubble grows 
up in size and strength. Later on, the upstream breakdown 
bubble also grows up in size and strength. The upstream 
breakdown bubble becomes larger and stronger than the 
downstream bubble, and the downstream bubble is con- 
vected through the duct exit at t = 78. Next the upstream 
breakdown bubble is convected downstreams and new 
breakdown bubbles appear behind the inlet shock. In the 
time range of t = 84-95, the mechanism of the evolution, 
merging, convection and shedding which is similar to that 
in the time range of t = 24-35 is repeated. In the time 
range of t = 96-120, the flow is similar to that of the time 
range of t = 37-78. At t = 123 and beyond, the whole 
process of vortex-breakdown-bubbles evolution, merging, 
convection and shedding is repeated. It is seen that the 
snapshots of streamlines at t = 124 and 130 are exactly 
similar to those at t = 3 and 17. Therefore, it is concluded 
that the vortex-breakdown mechanism for this flow case is 
periodic with a long period of time. Within this long cy- 
cle, short periodic cycles of vortex-breakdown develops. 
In summary, this flow case shows that several periodic 
modes of vortex breakdown develop, which correspond 
to different frequencies. 

Effect of Exit Boundary Conditions 

The Mach contours of the previous cases show that 
varying small portions of the duct exit become subsonic. 
Hence, extrapolating the components of the flowfield vec- 
tor from the interior cell centers for the duct exit boundary 
conditions is mathematically improper when it is viewed 
through the behavior of the characteristics at subsonic 
points. However, such boundary conditions could repre- 
sent typical physical conditions. Therefore, it is decided 
to examine the effect of using the Riemann-invariant-type 
boundary conditions at the subsonic points. This requires 
that four variables are extrapolated from the interior cell 
centers while the fifth must be specified at the exit We 
chose to specify the pressure at the duct exit at the sub- 
sonic points. 

For the flow case of Re = 100,000, we specified 
the exit pressure P* = P,* at the subsonic points. This 
type of boundary condition is enforced at t = 45 of the 
previous flow case (the solution of previous case at t 
s 45 serves as initial condition for the present case). 
Thereafter, the time-accurate integration with this type 
of boundary condition is continued. Figure 10 shows 
snapshots of the streamlines and Mach contours of this 
case. Although the solutions of this case in the time range 




of t = 46-61 is very similar to those of the previous case, 
a completely different flow for this case develops as of 
t ss 63. The two breakdown bubbles which are seen at 
t = 61 start moving in the downstream direction. They 
are continuously convected without any new generation 
of breakdown bubbles until they are completely shed 
through the small subsonic portion of the duct exiL In 
the meantime, the inlet shock disappears and the flow 
becomes supersonic throughout the duct 

Figures 11 and 12 show the results for Re = 20,000 
with P b = and P b = 2P«;, respectively. At t = 5, 
the breakdown bubble is exactly the same for these two 
cases, as it is seen from the streamline figures and the 
Mach- contour figures as well. At this instant of time, 
these solutions are the same as those of Fig. 6. It is 
obvious that the exit boundary conditions has not yet 
affected the upstream breakdown bubble. At the advanced 
time instants t = 18-20, we can see substantial differences 
between the solutions at each time instant of Figures 11, 
12 and 6. 

Thus, it is concluded that the exit boundary condi- 
tions have a substantial effect on the mechanism of break- 
down bubbles evolution, merging, convection and shed- 
ding. In particular, the solutions which are based on the 
Riemann-invariant-type boundary conditions are depen- 
dent upon the specified value of the pressure. It should 
be noted here that the Riemann-invariant-type boundary 
conditions, as it is well known is one dimensional and is 
based on inviscid analysis. Although the duct-exit flow 
is neither one-dimensional nor inviscid, the Riemann- 
invariant-type-boundary conditions are the best available 
non-reflective boundary conditions. Therefore, the ques- 
tion remains: what is the proper type of exit boundary 
condition that one needs to apply? Refer to Krause 15 for 
discussion on alternative boundary conditions. 

Effect of The Inlet Swirl Ratio 

In this flow case, the Reynolds number is kept at 
100,000 and the inlet swirl ratio is increased to 0.38. 
Figure 13 shows snapshots of the streamlines and Mach 
contours up to t = 30. Comparing the results of the present 
flow case with those of fi ■ 0.32, Figs. 8 and 9, we 
notice that the breakdown bubbles are larger and stronger. 
Moreover, in the time range of t = 17-30, the breakdown 
bubbles are oscillating around a mean position and a 
process of bubble generation, convection and shedding is 
taking place around a large breakdown bubble. The inlet 
shock has a very small amplitude oscillation around a 
mean position. This is a completely different mechanism 
from that of the case of Figs. 8 and 9. The computations 
of this flow case have not yet been carried out further 
in time. 


Concluding Remarks 

The unsteady, compressible NS equations are used 
for extensive computational study of supersonic quasi- 
axisymmetric vertex breakdown in a configured circular 
duct The quasi-axisymmetric solutions are time accu- 
rate. Several issues have been addressed in this paper. 
First, we have shown the effect of Reynolds number on 
the evolution, merging convection and shedding of vortex 
breakdown bubbles. Several modes of vortex breakdown 
have been obtained; a transient single-bubble breakdown, 
a transient multi-bubble breakdown and an unsteady peri- 
odic multi-bubble breakdown. These solutions have been 
obtained by using extrapolated flow conditions from the 
interior cell centers at the duct exit. Next, selected flow 
cases have been recomputed using a Riemann-invariant- 
type boundary conditions at the subsonic points of the 
duct exit It has been shown that substantial different 
solutions have been obtained and the question of what 
is the proper type of exit boundary conditions remains 
to be answered. Finally, the effect of the inlet swirl ra- 
tio has been investigated for the high Reynolds-number 
flow. It has been shown that a completely different mech- 
anism of vortex breakdown develops at advanced time 
instants. Work is underway to understand this mode of 
vortex breakdown. 

Acknowledgment 

For the first two authors, this work has been supported 
by the NASA Langley Research Center under Grant No. 
NAG-1-994. 

References 

1. Harvey, J. K., “Some Observations of The Vortex 
Breakdown Phenomenon,” Journal of Fluid Mech., 
Vol. 14, 1962, pp. 595-592. 

2. Sarpkaya, T., “Vortex Breakdown in Swirling Conical 
Flows,” AIAA Journal, Vol. 9, No. 9, Sept. 1971, 
pp. 1792^1799. 

3. Sarpkaya, T. f “Effect of The Adverse Pressure Gradi- 
ent on Vortex Breakdown," AIAA Journal, Vol. 12, 
No. 12, Dec. 1974, pp. 602-607. 

4. Leibovich, S., “Vortex Stability and Breakdown: Sur- 
vey and Extension ” AIAA Journal, Vol. 22, No. 9, 
Sept 1984, pp. 1192-1206. 

5. Escudier, M. P. and Zender, N. t “Vortex Flow 
Regimes" Journal of Fluid Mechanics, Vol. 115, 
1982, pp. 105-122. 

6. Grabowsld, W. J. and Berger, S. A., "Solutions of 
the Navier-Stokes Equations for Vortex Breakdown,” 
Journal of Fluid Mechanics, Vol. 75, Part 3, 1976, 
pp. 525-544. 

7. Hafez, M., Kuruvila, G. and Salas, M. D., "Numeri- 
cal Study of Vortex Breakdown,” Journal of Applied 
Numerical Mathematics, No. 2, 1987, pp. 291-302. 




8. Salas, M. D. and Kuruvila, G., "Vortex Breakdown 
Simulation: A Circumspect Study of The Steady, 
Laminar, Axisymmetric Model,” Computers and Flu- 
ids. Vol. 17, No. 1. 1989, pp. 247-262. 

9. Menne, S., “Vortex Breakdown in an Axisymmetric 
Flow,” AIAA 88-0506, January 1988. 

10. Wu, J. C. and Hwang, S., “Computational Study of 
Vortex Breakdown in Circular Tube,” AIAA 91-1820, 
June 1991. 

11. Menne, S. and Liu, C. H., “Numerical Simulation of a 
Three-Dimensional Vortex Breakdown,” Z. Flugwiss. 
Weltraumforsch. 14, 1990, pp. 301-308. 

12. Spall, R. E„ Gatski, T. B. and Ash, R. L„ “The Struc- 
ture and Dynamics of Bubble-Type Vortex Break- 
down,” Proc. R. Soc., London, A429, 1990, pp. 613- 
637. 

13. Breuer, M. and Hanel, D„ “Solution of The 3-D In- 
compressible Navier-Stokes Equations for the Simu- 
lation of Vortex Breakdown,” Eight GAMM Confer- 
ence, Delft, Netherlands, September 27-29, 1989. 

14. Krause, E„ “Vortex Breakdown: Physical Issues and 
Computational Simulation,” Third International Con- 
gress of Fluid Mechanics, Cairo, Egypt, January 1990, 
Vol. 1, pp. 335-344. 

15. Krause, E.,” The Solution of The Problem of Vortex 
Breakdown,” Invited paper, International Conference 
on Numerical Methods for Fluid Dynamics, Oxford, 
England, June 1990. 


16. Delery, J„ Horowitz, E., Leuchter, O. and Solignac, 
J. L., "Fundamental Studies of Vortex Flows,” La 
Recherche Aerospatiale, No. 1984-2, 1984, pp. 1- 
24. 

17. Metwally, O., Settles. G. and Horstman, C., “An Ex- 
perimental Study of Shock Wave/Vortex Interaction," 
AIAA 89-0082, January 1989. 

18. Cutler, A. D. and Levey, B. S., “Vortex Breakdown 
in a Supersonic Jet,” AIAA 91-1815, June 1991. 

19. Liu, C. H., Krause, E. and Menne, S., “Admissible 
Upstream Conditions for Slender Compressible Vor- 
tices,” AIAA 86-1093, 1986. 

20. Copening, G. and Anderson, J., “Numerical Solutions 
to Three-Dimensional Shock/Vortex Interaction at Hy- 
personic Speeds,” AIAA 89-0674, January 1989. 

21. Kandil, O. A. and Kandil, H. A., “Computation 
of Compressible Quasi-Axisymmetric Slender Vor- 
tex Flow and Breakdown,” IMACS 1st International 
Conference on Computational Physics, University of 
Colorado, Boulder, June 1990, pp. 46-51. Also in 
Computer Physics Communications, Vol. 65, North- 
Holland, March 1991, pp. 164-172. 

22. Meadows, K., Kumar, A. and Hussaini, M., “A Com- 
putational Study on the Interaction Between a Vortex 
and a Shock Wave,” AIAA 89-1043, April 1989. 

23. Kandil, O. A., Kandil, H. A. and Liu, C. H., "Compu- 
tation of Steady and Unsteady Compressible Quasi- 
Axisymmetric Vortex Flow and Breakdown,” AIAA 
91-0752, January, 1991. 




Figure 3. Streamlines and Mach contours for a swirling flow Figure 4. Streamlines for a swirling flow with transient 

with a transient breakdown, Moo = 1.75,/? = 0.32, multi-bubble breakdowns, = 1.75,/? = 0.32, 

R* = 4,000. Re = 10,000. 





























































































Figure 9. Mach contours for a swirling flow with unsteady 
multi-bubble breakdowns, Moo — 1.75, /J = 0.32, 
R* = 100,000. 





























Figure 13. Streamlines and Mach contours for a swirling flow with unsteady multi-bubble 
breakdowns, M* = 1.75, £ = 0.38, Re = 100,000. 
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Abstract 

The unsteady, compressible thin-layer and full 
Navier-Stokes equations are used to numerically simu- 
late steady and unsteady asymmetric, supersonic, locaily- 
conical flows around a 5°-semiapex angle circular cone. 
The main computational scheme used for the present com- 
putations is the implicit, upwind, flux-difference splitting, 
finite-volume scheme. Comparisons of the solutions us- 
ing the two sets of equations are presented for the flow 
asymmetry and its control. Computational studies are 
also presented to investigate the effects of the freestream 
Reynolds number and the locally-scaled Reynolds number 
on the flow asymmetry. These studies are carried out us- 
ing the full Navier-Stokes equations, Three-<limensional, 
asymmetric flow solutions are also presented for a 5°- 
semiapex angle cone of unit length and a cone-cylinder 
configuration. The three-dimensional solutions are ob- 
tained by using the thin-layer equations and short-duration 
transient side-slip disturbances along with a very fine grid. 

Introduction 

Highly swept, round and sharp-leading-edge wings 
and pointed slender bodies are common aerodynamic 
components to fighter aircraft and missiles. The study 
of vortex-dominated flow around these isolated aerody- 
namic components adds to our basic understanding of vor- 
tical flows under various conditions including unsteady 
vortex-dominated flows, vortex-shock interaction, asym- 
metric vortex flow and vortex breakdown. In this paper, 
we focus on the problem of asymmetric vortex flow about 
slender bodies in the high AOA range. The problem is 
of vital importance to the dynamic stability and control- 
lability of missiles and fighter aircraft 

The onset of flow asymmetry occurs when the relative 
incidence (ratio of angle of attack to nose semiapex angle) 
of pointed forebodies exceeds certain critical values. At 
these critical values of relative incidence, flow asymmetry 
develops due to natural and/or forced disturbances. The 
origin of natural disturbances may be a transient side slip, 
an acoustic disturbance, or likewise disturbance of short 
duration. The origin of forced disturbances is geometric 
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perturbations due to imperfections in the nose geometric 
symmetry or likewise disturbances of permanent nature. 
In addition to the relative incidence as one of the deter- 
minable parameters for the onset of flow asymmetry, the 
freestream Mach number, Reynolds number and shape of 
the body -cross sectional area are important determinable 
parameters. 

Kandil, Wong and Liu 1 used the unsteady, thin-layer 
Navier-Stokes equations along with two different implicit 
schemes to simulate asymmetric vortex flows around 
cones with different cross-section shapes. The numer- 
ical investigation was focused on a 5°-semiapex angle 
circular cone and local conical flow was assumed. The 
first computational scheme was an implicit, upwind, flux- 
difference splitting, finite-volume scheme and the sec- 
ond one was an implicit, central-difference, finite-volume 
scheme. Keeping the Mach number and Reynolds num- 
ber constants at 1.8 and 10 5 , respectively, the angle of 
attack was varied from 10° to 30°. At a = 10° , a steady 
symmetric solution was obtained and the results of the 
two schemes were in excellent agreement. At a = 20° 
and irrespective of the type or level of the disturbance, a 
steady asymmetric solution was obtained and the results 
of the two schemes were in excellent agreement. Two 
types of flow disturbances were used; a random round- 
off error or a random truncation-error disturbance and a 
controlled transient side-slip disturbance with short du- 
ration. For the controlled transient side-slip disturbance, 
the solution was unique, and for the uncontrolled random 
disturbance, the solution was also unique with the excep- 
tion of having the same asymmetry changing sides on the 
cone. At a = 30*, an unsteady asymmetric solution with 
vortex shedding was obtained, and the vortex shedding 
was perfectly periodic. Next, the angle of attack was 
kept fixed at 20* and the Mach number was increased 
from 1.8 to 3.0 with a step of 0.4. The solutions showed 
that the asymmetry become weaker as the Mach number 
is increased. The flow recovered its symmetry when the 
Mach number reached 3.0. Selected solutions of steady 
and unsteady asymmetric flows have also been presented 
for cones with elliptic and diamond cross-sectional areas. 
Passive control of the flow asymmetry has been tenta- 
tively demonstrated by using a fin on the leeward side of 
the body along the plane of geometric symmetry. 

Siclari* used the unsteady, Navier-Stokes equa- 
tions with a multi-grid, central-difference, finite-volume 
scheme to solve for steady asymmetric conical flows 
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around cones with elliptic, diamond and biparabolic sec- 
tions. He addressed steady-flow problems similar to those 
of the present authors in reference 1. He considered the 
flow around circular cones with semiapex angles of 5°, 
6°, T and 8° at an angle of attack of 20° and a Reynolds 
number of 10 5 . Varying the Mach number from 1.4 to 3.0 
with a step of 0.4, he showed that the flow recovered its 
symmetry as the Mach number increased. The higher the 
semiapex angle was, the lower the Mach number 
for the flow to recover it symmetry. Fixing the Mach 
number at 1.8, the angle of attack at 20°, the Reynolds 
number at I0 5 and the semiapex angle at 5°, he decreased 
the cross-section fineness ratio (ratio of width to length) 
for different cross-sectional shapes. He showed that the 
flow recovered its symmetry at a fineness ratio of 0.4 for 
the elliptic-section cone, at 0.6 for the biparabolic-section 
cone and at 0.6 for the diamond-section cone. 

In a very recent paper by Kandil, Wong, and Liu 3 , 
several issues related to the asymmetric flow solutions 
have been addressed. It has been shown that a unique 
asymmetric flow solution is obtained irrespective of the 
size of the minimum grid spacing at the solid boundary. 
The asymmetry could reverse sides due to the random 
nature of the disturbance. It has been also shown that 
for the same flow conditions and same section fineness 
ratio, diamond-section cones with sharp edges have less 
flow asymmetry than those of the elliptic-section cones. 
Moreover, it has been shown that passive control of flow 
asymmetry of diamond-section cones requires fin heights 
that are not necessarily equal to the local section width. 
On the other hand, passive control of flow asymmetry 
of circular and elliptic -section cones requires fins with 
heights that are, at least, equal to the local section width. 
Again, it was also shown that unsteady periodic asym- 
metric flow with vortex shedding has been predicted. 

In reference 4 by Kandil, Wong and Liu, several 
unsteady, asymmetric vortex flows with periodic vortex 
shedding for circular and noncircular section cones using 
the thin-layer Navier-Stokes equations w ere presented and 
studied. In reference 5 by Kandil, et al, the authors 
addressed the problem of asymmetric flow control using 
side strakes and the thin-layer Navier-Stokes equations. 
Two asymmetric flow cases have also been solved using 
the full Navier-Stokes equations. 

In the present paper, we address several issues related 
to the flow asymmetry around circular cones. Under 
the locally-conical flow assumption, steady and unsteady 
asymmetric solutions using the thin-layer and full Navier- 
Stokes equations are presented and compared. Also, 
control of flow asymmetry using vertical fins and side 
strakes are covered The solutions for the control cases 
are carried out using the thin-layer and full Navier-Stokes 
equations and their results are compared and discussed. 
Next, the effect of the frees tream Reynolds number on 
the flow asymmetry is studied under the locally-conical 
flow assumption using the full Navier-Stokes equations. 
This is achieved by keeping the axial station, at which the 


locally-conical solution is obtained, constant at the value 
of unity and changing the Reynolds number. The effect 
of scaled Reynolds number (Reynolds number based on 
the local axial station) on the flow asymmetry is also 
studied using the full Navier-Stokes equations. This is 
achieved by reducing the axial station and accordingly 
scaling the Reynolds number, the cross-section diameter, 
the grid fineness and the computational domain. Finally, 
three-dimensional flow asymmetry is investigated around 
a circular cone and a cone-cylinder configuration using 
the thin-layer Navier-Stokes equations. The purpose of 
this study is obtain flow asymmetry due to short-duration 
disturbances and to investigate the effect of Reynolds 
number, angle of attack and the cylinder afterbody on 
the flow asymmetry. 


Formulation 


Full and Thin-Layer Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady, 
compressible, full Navier-Stokes equations in terms of 
time-independent, body-conformed coordinates £ l , £ 2 and 
£ 3 is given by 


dQ dE m 

dt ^ d$ m 


3(Ev), 


= 0; m = 1 — 3, s = 


1-3 (1) 
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E ra = inviscid flux 
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(E V ) S = viscous and heat-conduction flux in 
direction 
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U m = <9^ ra u k (6) 

The first element of the three momentum elements of Eq. 
(5) is given by 
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The second and third elements of the momentum elements 
are obtained by replacing the subscript l, everywhere in 
Eq. (7), with 2 and 3, respectively. The last element of 
Eq. (5) is given by 


dk^HrVp - q k ) = 
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; p = 1 - 3 ( 8 ) 


The single thin-layer approximations of the full Navier- 
Stokes equations demand that we only keep the deriva- 
tives in the normal direction to the body, £ 2 , in the viscous 
and heat flux terms in Eqs. (1), (7) and ( 8 ). Thus, we let 

s = 2 for the term in Eq. (1) and s = 2 and n = 2 
in Eqs. (7) and ( 8 ). These equations reduce to 
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where 


= « 2 ^ 1 , * = l ^Sjf.w = a p ^ 2 u p (i2) 

The reference parameters for the dimensionless form 
of the equations are L, a^.L/a <x,P<x and for the 
length, velocity, time, density and molecular viscosity, 
respectively. The Reynolds number is defined as Re = 
PocVocL/^oc, where L is the body length, and the pres- 
sure, p, is related to the total energy per unit mass and 
density by the gas equation 


P = (7 “ 1 )P 
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The viscosity is calculated from the Sutherland law 

/i = T 3 / 2 ^ii^),C = 0.4317 (14) 

and the Prandtl number P r = 0.72. In Eqs. (1M12), the 
indicial notation is used for convenience. 


Boundary and Initial Conditions 

Boundary conditions are explicitly implemented. 
They include inflow-outflow conditions and solid- 
boundary conditions. At the plane of geometric symme- 
try periodic conditions are used for symmetric or asym- 
metric flow applications on the whole computational do- 
main (right and left domains). At the farfield inflow 
boundaries freestream conditions are specified since we 
are dealing with supersonic flows, while at the far-field 
outflow boundaries first-order extrapolation from the in- 
terior points is used. On the solid boundary, the no-slip 
and no-penetration conditions are enforced; ui = U 2 = U 3 
= 0 and the normal pressure gradient is set equal to zero. 
For the temperature, the adiabatic boundary condition is 
enforced on the solid boundary. 

For the passive control applications using vertical 
fins, double thin-layer Navier-Stokes equations are used, 
where one thin-layer is used normal to the body and 
another thin- lay er is used normal to the fin surface. For 
these applications, solid-boundary conditions are enforced 
on both sides of the fin. 

The initial conditions correspond to the uniform flow 
with m = U 2 = U 3 = 0 on the solid boundary. 


Computational Scheme 

The main computational scheme used to solve the 
thin-layer and full Navier-Stokes equations is an implicit, 
upwind, flux-difference splitting, finite-volume scheme. 
It employs the flux-difference splitting scheme of Roe. 
The Jacobian matrices of the in viscid fluxes, A $ = 
s = 1-3, are split into backward and forward fluxes ac- 
cording to the signs of the eigenvalues of the inviscid 
Jacobian matrices. Flux limiters are used to eliminate os- 
cillations in the shock region. The viscous and heat-flux 
terms are centrally differenced. The resulting difference 
equation is solved using approximate factorization in the 
(f 1 , £ 2 and 4 3 directions. The resulting computer program 
can be used to solve the thin-layer Navier-Stokes equa- 
tions and the full Navier-Stokes equations. This code 
is a modified version of the CFL3D which is currently 
called “FTNS3D” In this code, the implicit, flux- vector 
splitting, finite-volume scheme, which is based on the 
Van-Leer scheme 6 , can also be used instead of the flux- 
difference splitting scheme. 


Computational Applications and Discussions 

L Locally-Conical Asymmetric Flow Applications 

Locally-conical solutions of the thin-layer or full 
Navier-Stokes equations are obtained using one of two 
methods. In the first method, the governing equations 
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are transformed using the conical-coordinate transforma- 
tion. Invoking the conical flow conditions which re- 
quire that the flow variables be independent of the ra- 
dial distance (or axial distance, depending on the trans- 
formation) from the cone vertex, equating the radial dis- 
tance (or axial distance) which appears in the transformed 
equations to a constant (equals to unity in most of the 
present locally-conical solutions), the resulting equations 
are solved on one spherical (or cross-flow) surface. In the 
second method, the three-dimensional flow equations are 
solved on two spherical (or cross-flow) surfaces which 
are located in the very near proximity of a constant radial 
(or axial) distance. During the psuedo-time or accurate- 
time stepping, the flowfield vector is forced to be equal at 
the corresponding grid centers on the two surfaces. This 
method is used in the present paper to obtain locally- 
conical solutions. The resulting solutions from these two 
methods are the same locally-conical solutions. These 
solutions correspond to the specified radial (or axial) dis- 
tance and hence they change as we change the radial (or 
axial) distance. The reason behind that is simply because 
the transformed equations, according to the first method, 
are not self-similar and hence they are not globally coni- 
cal. This is shown below by developing the transformed 
equations of the first method. Considering the unsteady, 
compressible, Navier-Stokes equations in the Cartesian 
coordinates. 


introducing the conical coordinates, 
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and using the chain rule, Eq. (15) is transformed to 
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The conical flow condition requires the flow variables be 
independent of the coordinate m (radial distance). Invok- 
ing this condition in Eq. (17) by dropping the derivatives 
with respect to m< Eq. (17) reduces to 

slr + 4; (£ - 6 ' ) ' + ^ (f - £ ' ) - +2(i - i - )=0 

(19) 

It is obvious that the unsteady term includes 773. More- 
over, the viscous terms 4^, and I, include m> and 


hence Eq. (19) is not self-similar. The explicit depen- 
dence of the viscous terms on rj, can be shown through 
one of the elements of these vectors. For example, we 
consider 
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Thus, the unsteady term and viscous terms are scaled by 
the radial distance 773 and Eq. (19) does not represent a 
globally-conical flow. The best to be done to make use 
of this equation is to select a constant value for 773, and 
solve the resulting equation for what we call “locally- 
conical flow”. If 73 is assigned another constant value, the 
resulting equation will have another scale for the unsteady 
term and viscous terms. It is concluded that Eq. (19) 
becomes globally conical if the unsteady term and viscous 
terms vanish, and hence only the steady Euler equations 
are globally conical. 

Next, we present comparison of the solutions using 
the thin-layer and full Navier-Stokes equations. We also 
present steady flow solutions using different Reynolds 
numbers keeping 773 = 1 and steady solutions using 773 
= 0.5, 0.25, 0. 1 and Re = 50,000; 25,000; 10,000; re- 
spectively. 


Steady Asymmetric Flow (a = 20°, M* = 1.8, R* = 10 5 ) 

Supersonic flow around a 5° -semiapex angle circular 
cone at 20° angle of attack, freestream Mach number 
of 1.8 and Reynolds number of 10 5 is considered. The 
thin-layer and full Navier-Stokes equations are used to 
solve for asymmetric flow on a grid of 161 x81 points in 
the circumferential and normal directions, respectively. 
The grid is generated by using a modified Joukowski 
transformation with a geometric series for grid clustering 
in the normal direction. The minimum grid spacing in 
the normal direction at the solid boundary is 10" 4 . The 
maximum radius of the computational domain is 21 r, 
where r is the radius of the circular cone at the axial 
station of unity. Figure 1 shows the results of thin- 
layer solution. The color graphics figure shows six snap 
shots for the evolution of the steady asymmetric solution 
during 10,000 iteration steps. The number from 1-6 on 
the total-pressure-loss contours correspond to the number 
from 1-6 on the logarithmic residual-iterations curve. It 
is seen that the solution is symmetric during the first 3000 
iteration when the residual drops to machine zero (No. 1). 
Then, the solution becomes slightly asymmetric as the 
residual grows up (Nos. 2, 3). As the residual drops 
again to machine zero, the flow asymmetry becomes 
strong until it reaches to the stable asymmetric solution 
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(Nos. 4, 5, 6). The final total -pressure -loss contour and 
the corresponding surface-pressure coefficient are also 
shown. 

The full Navier-Stokes solution is shown in Fig. 2, 
The residual-error drops 4.5 orders of magnitude in 2,000 
steps, increases one order of magnitude after a total of 
3,000 steps and then drops to machine zero in a total 
number of iteration steps of 6,000. It is clear that the 
full Navier-Stokes equations produce the asymmetric so- 
lution faster than the thin-layer Navier-Stokes equations. 
The total-pressure-loss contours show that the full Navier- 
Stokes solution produces thicker shear layers than those 
of the thin-layer solution. More contour resolution in the 
vortex cores is produced by the full Navier-Stokes so- 
lutions than that of the thin-layer solution. Finally, the 
free-shear layer on the body right-side of the full Navier- 
Stokes solution is shorter than that of the thin-layer so- 
lution. However, the Cp figures of the two solutions are 
exactly the same. 

Asymmetric Flow Control Using Vertical 
Fins (h = 0.5r, h = r) 

Figures 3 and 4 show comparisons of the thin-layer 
and full Navier-Stokes solutions for the control of the 
flow of the preceding case. Two vertical fins of heights h 
= 0.5r and r are placed in the leeward plane of geometric 
symmetry, where r is the cone local radius. Figures 3 and 
4 show the total-pressure-loss contours and the surface- 
pressure coefficient. The thin-layer solutions are shown 
on the left and the full Navier-Stokes solutions are shown 
on the right. With h = 0.5r and using the thin-layer 
equations, two vortex cores which are connected to each 
other and to the body through free-shear layers, develop 
from the left side of the body. From the right side of 
the body, a free-shear layer develops and crosses over 
the fin to the left side of the body. It produces two 
vortex cores; one at each comer of the body-fin juncture 
with secondary separation below them. This case has 
been solved accurately in time but it does not show any 
vortex shedding or unsteadiness. When the fin height is 
increased to h = r, perfect flow symmetry is obtained. 

With h = 0.5r and using the full Navier-Stokes equa- 
tions, a mirror image of the vortex system and flow sepa- 
rations of the thin-layer solutions is obtained. This is due 
to the random nature of the disturbance. However, the 
freeshear layer and vortex cores on the right are thicker 
and shorter in height than those of the thin-layer solution. 
However, the C p coefficients of the two sets of equations 
are exact mirror image of each other. With h = r, the full 
Navier-Stokes equations produce the exact same solution 
as that of the thin-layer equations. The reason behind the 
flow asymmetry with h = 0.5r is that the free-shear layer 
from the right-hand side of the body is still higher than 
the fin height, which allows the flow disturbances from 
the right and left side to interact 


Asymmetric Flow Control Using Side 
Strakes (h = OJr) 

Figure 5 shows a comparison of the thin-layer and 
full Navier-Stokes solutions for the control of the flow 
of the case of Fig. 1. Side strakes of height h = 0.3r, 
when r is the local radius of the cone, are used for this 
purpose. The grid used is generated using a hyperbolic 
grid generator with transfinite grid interpolation to refine 
the grid at the strake sharp edge. While the thin-layer 
solution shows perfect symmetric flow, the full Navier- 
Stokes solution shows slightly asymmetric flow. This 
can be seen by comparing both the total-pressure-loss 
contours and the surface-pressure coefficient of the two 
solutions. Both solutions are obtained using a grid of 
161x81 with the same transfinite interpolation. Both 
surface-pressure coefficients show a jump in the pressure 
value at the leading edges of the strakes. It should be 
noted that side-strakes produce higher lift than that of the 
vertical fin control. 

Unsteady Asymmetric Flow (a = 30°, 

= L8, R. = 10 5 ) 

Figures 6 and 7 show comparisons of the full Navier- 
Stokes solutions and the thin-layer solutions for unsteady 
asymmetric flow with vortex shedding. The solutions 
are obtained by using the full Navier-Stokes equations 
with a grid of 161x81 and the flux difference splitting 
(FDS) scheme, thin-layer equations with a grid of 16 1 x 8 1 
and the FDS scheme, thin-layer equations with a grid 
of 161x81 and the flux-vector splitting (FVS) scheme, 
and full Navier-Stokes equations with a grid of 241x81 
and the FDS scheme. Typical of all the four solutions, 
the residual error curves show a drop in the error fol- 
lowed by a transient response and ending by a periodic 
response. The lift coefficient curves show the final peri- 
odic response. The first three solutions exactly show the 
same number of time steps (At = 1CF 3 ) of 700 for one-half 
of the cycle of vortex shedding. The full Navier-Stokes 
solution on the fine grid (Fig. 7) shows a slightly bigger 
number of time steps (At = 1CT 3 ) of 725 for one-half of 
the cycle of vortex shedding. This is expected of the full 
Navier-Stokes equations on the fine grid since the flow 
has more real damping and viscous-force resolution than 
the other three solutions. In the fourth case, the damping 
will elongate the period of shedding and hence it reduces 
the frequency of shedding. The frequency of shedding 
of the first three solutions is 4.488 and the frequency of 
shedding of the fourth solution is 4.333. 

Asymmetric Flow Control Using a Vertical Fin (h = r) 

Figure 8 shows the full Navier-Stokes solution of 
asymmetric flow control using a vertical fin of h = r. This 
is the flow considered in the preceding case. It is seen 
that the vertical fin height is not enough to eliminate the 
flow asymmetry. The flow is still asymmetric but steady. 
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Effect of Decreasing the Reynolds Number Keeping 
the Axial Distance Fixed (x = 1, a = 20°, 

= 1.8, R € = 20,000; 15,000; 10,000) 

Figure 9 shows the results of a study to investigate 
the effect of reducing the Reynolds number on the flow 
asymmetry. The full Navier-Stokes equations are used on 
a grid of 241 x 81. As the Reynolds number is decreased 
keeping x = 1, the flow asymmetry decreases but it does 
not completely diminish. In this study, the computational 
domain is kept fixed in size and the minimum grid spacing 
at the solid boundary is kept at lO -4 . It should be stressed 
here that with decreasing the Reynolds number to the 
values considered, only the full Navier-Stokes solutions 
on a relatively fine grid will be valid. Strictly speaking, 
the cut-off Reynolds number for flow asymmetry is below 
10 , 000 . 

Effect of Decreasing the Reynolds Number and 
the Axial Distance (x = 0.5, 0.25, 0.1, a = 20°, 

Moc = 1.8, R* = 50,000; 25,000; 10,000) 

Figure 10 shows the results of a study to investigate 
the effect of reducing the axial distance and proportionally 
the Reynolds number, the computational domain and the 
minimum grid spacing at the solid boundary. In other 
words, we are investigating whether the flow is self- 
similar or not. The solutions shown in the figure ensure 
that the flow is not self similar as the problem is scaled 
at each axial station. However, in another numerical 
study, we have shown that if the Reynolds number is 
increased by the same ratio as that of decreasing the 
axial station (e.g. if x = 0.5 then R* = 200,000 such 
that xR e = 100,000) and the problem is solved for the 
corresponding reduced radius of the cone section and 
its reduced computational domain, the asymmetric flow 
solution is self similar at any section. 

Unsteady Asymmetric Flow (a = 35°, 

Moc = U, Re = 10 5 

Figure 11 shows snap shots of the unsteady asym- 
metric flow of the full Navier-Stokes equations for the 
same circular cone as the angle of attack is increased 
to 35°. The residual error curve and the lift coefficient 
curve show the same typical responses as those of Figs. 6 
and 7. The period of shedding is smaller and equals to 
550xl0~ 3 = 0.55 and the corresponding shedding fre- 
quency is 11.424. 

n. Three-Dimensional Flow Applications 

In Figs. 12—19, we present the results of three- 
dimensional asymmetric flow solutions of the thin-layer 
equations on a grid of 161 x81 x65 in the circular, normal 
and axial directions, respectively. There are several issues 
to be answered through the present computational study. 
First, for the same circular cone and for the same flow 
conditions (angle of attack, Mach number and Reynolds 
number), will the three-dimensional flow solution be the 
same as that of the locally-conical solution of Fig. 1? If 


the answer is negative, the next question to address is: Is 
there a length scale which relates the three-dimensional 
solution to the locally-conical solution? The second issue 
to be addressed is the effect of the Reynolds number, the 
angle of attack and the cyhndrical-afterbody length on the 
three-dimensional flow asymmetry? 

To address these issues, a 5°-semiapex angle circu- 
lar cone of unit length (cone length is the characteris- 
tic length) is considered. The three-dimensional grid of 
161x81x65 is generated by using a modified Joukowski 
transformation at axial stations. The grid is clustered in 
the normal direction of the body using a geometric series 
with the minimum grid spacing of l(faat the vertex and 
10“*" at the axial station of unit length. A typical grid is 
shown in Fig. 12. With the flow conditions set at a = 
20°, Moc = T8 and R« = 10 5 , which are the same con- 
ditions as those of the locally conical flow of Fig. 1, the 
three-dimensional solution produced a symmetric flow, 
unlike the local-conical solution which produces asym- 
metric steady flow. The reason for the difference is weU 
understood since the local-conical solution is obtained at 
an axial station of x = 1.0, Hence, a length scale is 
involved in the Reynolds number, as can be seen from 
the analytical conical equation for steady viscous flow, 

9 Eq* (19). Next, the search is directed at obtaining asym- 
metric flow solutions for the three-dimensional core flow. 
In Fig. 13, we show the solution for the same cone at a 
= 40°, M<x » 1.4 and R« - 4x 10 6 . It is seen that the 
solution is asymmetric and is nearly self-similar over a 
long axial distance of the cone length. 

Next, the Reynolds number is increased to 6x 10 6 and 
8x10* keeping the other flow conditions constant at a 
= 40° and Moc = 1.4. Figures 14 and 15 show strong 
asymmetric-flow solutions with already shed vortices. It 
should be noticed that the flow asymmetry changes sides 
as we move in the downstream direction. Hence, we have 
spatial asymmetric vortex shedding which is qualitatively 
similar to the temporal asymmetric vortex shedding of the 
locally-conical flow solutions of Figs. 6, 7 and 11. The 
flow instability is of spatial type. 

In an attempt to address the issue of the effect of the 
cylindrical afterbody on the flow asymmetry, a cylindrical 
afterbody of unit length is added to the unit conical 
forebody. The flow configuration is solved for the flow 
conditions of the isolated cone of Fig. 13. The results 
are shown in Fig. 16. Comparing the total-pressure- 
loss contours of Figs. 13 and 16, we see that the flow 
asymmetry is stronger for the cone-cylinder configuration 
in comparison with that of the cone alone. It should be 
noted that inside the shock cone surrounding the cone- 
cylinder configuration, subsonic flow exists and hence the 
downstream boundary has an upstream effect Figures 17 
and 18 show the total-pressure-loss contours and surface- 
pressure coefficients in cross-flow planes for the cases of 
Figs. 15 and 16. 
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Finally, the angle of attack is increased to SO 4 keep- 
ing the Mach number and Reynolds number constants at 
1.4 and 8x 10 6 , respectively. The solution is obtained ac- 
curately in time. In Fig. 19, we show two snap shots of 
the solution at the time steps of n = 10,116 and 11,818 
which correspond to t = 0.10116 and 0.11818, respec- 
tively, since At = 10' 5 . Comparing the two solutions, one 
can see vortex shedding. This case of three-dimensional 
unsteady vortex shedding is computationally very expen- 
sive since the stable time step is 10" 5 , and hence it was 
not completed. 

Concluding Remarks 

In the present paper, the unsteady, compressible thin- 
layer and full Navier-Stokes equations have been used 
to solve for steady and unsteady asymmetric flows and 
their passive control around a 5°-semiapex angle cone. 
For the steady locally-conical asymmetric solutions, we 
have shown that the full Navier-Stokes solutions pro- 
duce thicker and shorter ffee-shear layers than those of 
the thin-layer solutions. For the unsteady locally-conical 
asymmetric solutions, we have shown that the full Navier- 
Stokes solutions on a fine grid produce longer periods of 
vortex shedding and hence smaller shedding frequencies 
than those of the thin-layer solutions. Next, we addressed 
the effects of reducing the Reynolds number on the flow 
asymmetry. First, the Reynolds number is reduced keep- 
ing the axial distance constant. We have shown that the 
cut-off Reynolds number for the flow asymmetry is less 
than 10,000. Second, the Reynolds number is reduced 
along with the proportional reduction of the axial distance, 
the computational domain and the minimum grid spacing. 
We have shown that the flow asymmetry disappears be- 
low the Reynolds number value of 25,000. Moreover, 
we have shown that the flow is not globally conical. Fi- 
nally, the three-dimensional flow problem is addressed for 
the same cone and a cone-cylinder configuration. Flow 
asymmetry has been obtained using short-duration dis- 
turbances. The flow asymmetry becomes stronger as the 
Reynolds number and the angle of attack are increased. 
It also becomes strong due to the addition of a cylindrical 
afterbody. Unsteady flow asymmetry has also been ob- 
tained. It has also been noticed that for certain flow condi- 
tions, the flow asymmetry shows spatial vortex shedding 


which is qualitatively similar to the temporal vortex shed- 
ding of the unsteady locally-conical asymmetric flow. 
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thin-layer equations. 
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Figure 8. Passive flow control using a vertical fin, h = r, full Navier-Stokes solution. 
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Figure 9. Effect of decreasing the Reynolds number keeping the axial station constant, x = l, 241 x 81, 
full Navier-Stokes solutions. 
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Figure 10. Effect of decreasing the Reynolds number and the axial station, R« = 50,000; 25,000; 10,000; 


x = 0.5, 0.25, 0.1; 241x161, full Navier-Stokes solutions. 
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Abstract 

The unsteady, compressible, thin-layer Navier-Stokes 
equations are used to obtain three-dimensional, asym- 
metric, vortex-flow solutions around cones and cone- 
cylinder configurations. The equations are solved us- 
ing an implicit, upwind, flux-difference splitting, finite- 
volume scheme. The computational applications cover 
asymmetric flows around a 5°semi-apex angle cone of 
unit length at various Reynolds number. Next, a cylin- 
drical afterbody of various length is added to the conical 
forebody to study the effect of the length of cylindrical 
afterbody on the flow asymmetry. One of the compu- 
tational solutions has been validated by comparing the 
computed surface pressure with those of the experimen- 
tal data. All the asymmetric flow solutions have been 
obtained by using a short-duration side-slip disturbance. 

Introduction 

In this paper, the problem of asymmetric vortex-flow 
around three-dimensional cone and cone -cylinder body is 
addressed. This problem has received considerable at- 
tention by researchers in the computational fluid dynam- 
ics area 1 ' 7 and by researchers in the experimental fluid 
dynamics area 6-13 . The problem is of vital importance 
to the dynamic stability and controllability of missiles 
and fighter aircraft When flow asymmetry develops, it 
produces side forces, asymmetric lifting forces and cor- 
responding yawing, rolling and pitching moments that 
might be larger than those available by the control sys- 
tem of the vehicle. The onset of flow asymmetry oc- 
curs when the relative incidence (ratio of angle of attack 
to nose semi-apex angle) of pointed forebodies exceeds 
certain critical values. At these critical values of rela- 
tive incidence, flow asymmetry develops due to natural 
and/or forced disturbances. The origin of natural dis- 
turbances may be a transient side-slip, an acoustic dis- 
turbance, or similar disturbance of short duration. The 
origin of forced disturbances is geometric perturbations 
due to imperfections in the nose geometric symmetry or 
similar disturbances of permanent nature. In addition to 
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the relative incidence as one of the influential parameters 
for the onset of flow asymmetry, the frees tream Mach 
number, Reynolds number and shape of the body-cross 
sectional area are also important parameters. 

In several recent papers by the present authors 1 * 4 , the 
unsteady, thin-layer, compressible Navier-Stokes equa- 
tions have been used to simulate steady and unsteady, 
asymmetric vortex flows, including their passive control, 
around cones with different cross-sectional shapes. The 
emphasis of these papers was extensive computational 
studies of the parameters which influence the asymmet- 
ric flow phenomenon and its passive control. Since the 
computational cost associated with the solution of three- 
dimensional-flow problems with reasonable flow resolu- 
tion is very expensive, all the computational solutions 
were obtained using a Iocally-conical flow assumption. 
Such an assumption reduces the problem solution to that 
on two conical planes, which are in close proximity of 
each other, and hence it reduces the computational cost 
by an order of magnitude. Moreover, such solutions still 
provide extensive understanding of the flow physics since 
one can use very fine grids for reasonable flow resolution. 
These studies showed that asymmetric flow solutions 
were unique irrespective of the type of flow disturbance; 
a random disturbance in the form of a machine round-off 
error or a controlled disturbance in the form of a short- 
duration side-slip disturbance. Unsteady asymmetric flow 
solutions with perfectly periodic vortex shedding were 
successfully simulated, and the solutions were unique ir- 
respective of the computational scheme used. We also 
showed that as the Mach number was increased, the flow 
asymmetry was decreased and as the Reynolds number 
was increased, the flow asymmetry was increased. The 
cross-sectional shape of the cone has been shown to be a 
very influential parameter on the flow asymmetry. Circu- 
lar sections produced very strong flow asymmetry and di- 
amond sections produced relatively- weaker flow asymme- 
try. Passive control of the flow asymmetry was demon- 
strated by using vertical fins of different heights along 
the leeward plane of geometric symmetry and by using 
thin and thick side strakes with different orientations. It 
was also shown that side-strakes control is more practical 
than the vertical-fin control since it was effective over 
a wide range of angle of attack and provided additional 
lifting force. In a later paper, by the present authors 3 , 
the full Navier-Stokes solutions were compared with the 
thin-layer Navier-Stokes solutions. It was shown that the 
full Navier-Stokes solutions produced thicker free-shear 
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layers and more vortex-core resolution as compared with 
those of the thin-layer Navier-Stokes equations. In refer- 
ence 5 , a few tentative three-dimensional flow solutions 
were also presented. 

In an attempt to simulate asymmetric vortex flow 
around an ogive-cylinder body at an angle of attack of 
40°, a freestream Mach number of 0.2 and a freestream 
Reynolds number of 200,000; Degani and Schiff 5 used 
the unsteady, thin-layer, Navier-Stokes equations along 
with an implicit scheme which is second-order accurate in 
time. The scheme uses central-differencing in the cross- 
flow plane and upwind flux-vector splitting in the stream- 
wise direction. By introducing a forced asymmetric dis- 
turbance near the body nose in the form of a small surface 
jet, asymmetric flow solution was obtained. When the jet 
was turned off, the flow asymmetry dissipated and the 
flow recovered its symmetry. 

In a later paper by Degani 7 , the same computational 
scheme was used to solve for the flow around the same 
ogive-cylinder body over a wide range of angle of attack; 
a = 20° - 80°. His numerical experiments focused on 
investigating the origin of vortex asymmetry. Based on 
his results, he suggested that the flowfield around slender 
bodies could be divided into three main groups depending 
on the angle of attack range. This range might change 
by ±10°, depending on the flow conditions. In the range 
0° < a < 30°, the flow was symmetric and introduction 
of small disturbances near the nose had a small effect on 
the flow symmetry. In the second range, 30° < a < 60°, 
the flow became steady asymmetric upon introduction of 
a space-fixed forced disturbance near the nose. The level 
of asymmetry was a function of the location and size of 
the forced disturbance, and for large size disturbances, the 
asymmetry became unsteady with very high frequency. 
However, when the disturbance was removed the flow 
recovered its symmetric shape. He attributed the origin 
of asymmetry to a convective-type-instability mechanism. 
In the very high range, 60° < a < 80°, the flow became un- 
steady with vortex shedding upon introduction of a small 
transient disturbance with short duration. He attributed 
the origin of flow unsteadiness and vortex shedding to 
an absolute-type-instability mechanism. In that range of 
angle of attack, he also showed that the convective-type- 
instability mechanism was possible upon introduction of 
a space-fixed disturbance near the nose. Although this 
investigation revealed good tentative conclusions, there 
are several remaining questions to be addressed, which 
are related to the scheme dissipative effects, particularly 
in the cross-flow planes, and the grid fineness and its 
resolution of the disturbance growth. 


slip. With guidance from the locally-conical solutions, the 
angle of attack is varied between 30° and 50°, the Mach 
number is varied between 1.4 and 1.8 and the Reynolds 
number is varied between 10 5 and 8x10 s searching for 
asymmetric flow solutions due to a short-duration distur- 
bance. Next, the flow conditions are fixed and a cylindri- 
cal afterbody is added to the same cone of unit length to 
study the effect of the length of the cylindrical afterbody. 
The computational results have also been verified using 
comparison of the surface-pressure coefficient with that 
of the experimental data of Landrum 8 . 


Formulation 


Thin-Layer Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady, 
compressible, thin-layer Navier-Stokes equations in terms 
of time -independent, body-conformed coordinates £ 2 
and £ 3 is given by 
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(E v ) j = viscous and heat— conduction flux in 
direction 

ak£ 2 (u„fti - qit)]*; k = 1 - 3, n = 1 - 3 (5) 

u m = a^ m u M (6) 

The first element of the three momentum elements of Eq. 
(5) is given by 


In the present paper, we focus on the three- 
dimensional asymmetric flow problem. In particular, 
we address several important issues concerning the flow 
asymmetry around three-dimensional bodies. First, the 
three-dimensional asymmetric flow around a 5°semi-apex 
angle cone of unit length is considered in response to a 
short-duration disturbance in the form of a transient side 


The second and third elements of the momentum elements 
are obtained by replacing the subscript 1, everywhere in 


2 



Eq. (7). with 2 and 3, respectively. The last element of 
Eq. (5) is given by 


dk£ 2 (u„TTcn - qic) = 


+4 
+ 


Re 
l_d 
2d? 

1 <?(a 2 ) 


v W 


! (UpUp) 


(7 - l)Pr <?£ 2 J 


; P = 1 - 3 (8) 


where 


4> = = W = a„^u n (9) 


In Eqs. (l)-(9), the dimensionless variables are referenced 
to their appropriate freestream values. The dimensionless 
density p, Cartesian velocity components ui, U2 and 113, 
total energy per unit mass, e, dynamic viscosity, p and 
speed of sound, a, are defined as the ratio of the corre- 
sponding physical quantities to those of the freestream; 
namely, p^, a^, p^ and a^; respectively. The 

pressure, p, is non-dimensionalized by and is re- 

lated to the total energy for a perfect gas by the equation 
of state 


P = (7 - 1 )p (e ~ ^ u i u i) ; j = 1 - 3 (10) 

where 7 is the ratio of specific heats and its value is 1.4. 
The viscosity, p> is calculated from the Sutherland’s law 

fi = T 3/2 (iil) , c ss 0.4317 (11) 

where T is the temperature which is non-dimensionalized 
by T«;. The Prandtl number, ?, is fixed at 0.72. The 
Reynolds number is defined as R« = p< x U 0C L/^, x and 
the characteristic length, L, is chosen as the length of the 
body. 

In Eqs. (1)-(10), the indicial notation is used for con- 
venience. The subscripts k, n, p and j are summation 
indices, the superscript or subscript s is a summation in- 
dex and the superscript or subscript m is a free index. 
The partial derivative ^ is referred to by 

Boundary and Initial Conditions 

Boundary conditions are explicitly implemented. 
They include the inflow-outflow conditions and the 
solid-boundary conditions. At the plane of geomet- 
ric symmetry, periodic conditions are used. Since the 
freestream Mach number is supersonic and the inflow- 
outflow boundaries are also supersonic, freestream condi- 
tions are specified at the inflow boundaries and first-order 
extrapolation of the flow variables is used at the outflow 
boundaries. The conical shock enclosing the body is cap- 
tured as part of the solution. On the solid boundary, the 
no-slip and no-penetration conditions are enforced; u< = 
U2 = uj = 0, and the normal pressure gradient is set equal 


to zero. For the temperature, the adiabatic boundary con- 
dition is enforced at the solid boundary. 

The initial conditions correspond to the freestream 
conditions with ui = u? = U3 = 0 on the solid boundary. 
The freestream. conditions are given by 


Poc — — T^c — 1, 

Ujoc = XL* cos a cos/?, 
U’cc — -M* sin ,3, 


uj* = M<x sin a cos /?, 


Poc = 1/7. 
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7(7-1) 



( 12 ) 


where a is the angle of attack and j3 the side slip angle. 


Computational Scheme 

The implicit, upwind, flux-difference splitting finite- 
volume scheme is used to solve the unsteady, compress- 
ible, thin-layer Navier-Stokes equations. The scheme uses 
the flux-difference splitting scheme of Roe which is based 
on the solution of the approximate Riemann problem. 
In the Roe scheme, the inviscid flux difference at the 
interface of computational cells is split into two parts; 
left and right flux differences. The splitting is accom- 
plished according to the signs of the eigenvalues of the 
Roe averaged- Jacobian matrix of the inviscid fluxes at 
the cell interface. The min-mod flux limiter is used to 
eliminate oscillations in the shock region. The viscous- 
and heat-flux terms are linearized and the cross-derivative 
terms are eliminated in the implicit operator. The vis- 
cous terms are differenced using a second-order accurate 
central differencing. The resulting difference equation is 
approximately factored and is solved in three sweeps in 
tire ? and ? directions. The computational scheme 
is coded in the computer program “CFL3D.” 


Computational Applications and Discussions 

In the present computational applications, we consider 
the three-dimensional solutions of the unsteady, com- 
pressible, thin-layer Navier-Stokes equations for asym- 
metric vortex flows around a circular cone and circu- 
lar cone -cylinder configurations. There are several is- 
sues to be addressed through the present study. First, 
for the same circular cone of 5°-semi-apex angle 1 and 
for the same flow conditions and source of disturbance, 
will the three-dimensional flow solution be the same as 
that of the locally-conical flow solution? If the answer 
is negative, the next question to address is; Is there a 
length scale which relates the three-dimensional solution 
to the locally-conical solution? The second issue to be 
addressed is the effect of Reynolds number on the flow 
asymmetry? The third issue to be addressed is the effect 
of the cylindrical-afterbody length on the flow asymme- 
try. Finally, we address the question of code and grid val- 
idation by comparing the results of an asymmetric flow 
solution with those of the experimental data 1 . 
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Circular Cone 

A 5°-semi-apex angle circular cone of unit length 
(cone length is the characteristic length) is considered to 
address the first two issues mentioned above. This is the 
same circular cone which was considered by the authors 
in Ref. 1 for the locally-conical flow solutions. A three- 
dimensional grid of 161 x 81 x65 in the wrap around, nor- 
mal and axial directions, respectively, is generated by us- 
ing a modified Joukowski transformation at axial stations. 
The grid is clustered algebraically in the normal direction 
of the body using a geometric series with minimum grid 
spacing of lO^ 5 at the cone vertex and 10“ 5 at the axial 
station of unit length. A typical grid is shown in Fig. 1. 
The cross-flow grid size of 161x81 is the same grid size 
which was used for the locally-conical flow solutions of 
Ref. 1. 

With the flow conditions set at a = 20°, M* = 1.8 and 
R c = 10 5 , which are the same conditions as those of the 
locally-conical flow of Ref. 1, the three-dimensional solu- 
tion produces a symmetric steady flow, unlike the locally- 
conical solution which produces asymmetric steady flow. 
The reason for the difference is well understood since 
the locally-conical solution is obtained at an axial sta- 
tion of x = 1.0. Hence, a length scale is involved in the 
Reynolds number, as can be seen from the analytical con- 
ical equation for steady viscous flow. Next* the search is 
directed at obtaining asymmetric flow solutions for the 
three-dimensional cone flow. In Fig. 2, we show the so- 
lution in the form of total-pressure loss for the same cone 
at a = 40°, Moc = 1.4 and R« = 4X10 6 . It is seen that 
the solution is asymmetric and is nearly self-similar over 
a long axial distance of the cone length. This solution 
is obtained using a short-duration side-slip disturbance. 
When the residual error drops four orders of magnitude, 
a side-slip disturbance of 0 = 2° is applied for 100 it- 
eration steps, then it is removed. Thereafter, the pseudo 
time stepping is continued until the residual error drops 
again four to five orders of magnitude and a stable asym- 
metric solution is obtained. It should be noted that a ma- 
chine round-off error type of random disturbance cannot 
be used to obtain three-dimensional asymmetric solutions 
since the residual error never drops to machine zero. 

Next, the Reynolds number is increased to SxlO 6 
and 6X10 6 keeping the other flow conditions constant 
at a = 40 and 1.4. Figures 3 and 4 show the total- 
pressure-loss solutions for these cases. Figure 3 shows 
that the asymmetry of the vortical flow gets strong and the 
self similarity of the flow asymmetry is substantially lost. 
However, it is noticed that the flow asymmetry does not 
change sides as the solution develops in the downstream 
direction. Figure 6, which corresponds to the R* of 
fixlO 6 , shows that the flow asymmetry changes siiies 
as the solution develops in the downstream direction. 
Moreover, it is noticed that shed vortices exist in the flow, 
A close study of the solutions between the shown fourth 
cross-flow plane and seventh cross-flow plane reveals 
that the flow asymmetry changes from the right side 


(fourth cross-flow plane) to the left side (seventh cross- 
flow plane). The solutions on these two planes are nearly 
scaled mirror-images of each other. The present spatial 
flow asymmetry is qualitatively similar to the temporal 
flow asymmetry of the locally-conical flow solution of 
Ref. 1 (see Figure 8 of the present paper). 

Figures 5 and 6 show front and rear side views of 
the limiting streamlines for the cases of R* = 4x 10 6 and 
R c = 6x 10 6 , respectively. By comparing the streamlines, 
lines of separation and attachment of the front and rear 
side views in each figure, it is noticed that the flow asym- 
metry exists on the boundary and becomes stronger as 
the Reynolds number increases. As the Reynolds number 
increases, the separation and reattachment lines changes 
from radial straight lines to highly curved lines. 

Figure 7 shows the total-pressure-loss solution for the 
same cone for a higher Reynolds number, R« = 8x 10 6 . 
The asymmetry of the vortex flow becomes much stronger 
as compared with the previous cases of Figs. 2 - 4 , By 
comparing the solution of this case with that of the R c 
= 6xl0 6 , it is noticed the flow asymmetry of the case 
with high R* changes sides along a shorter axial distance 
(third and fifth cross-flow planes) in comparison with that 
of the low R«. Moreover, the flow asymmetry of the 
case with high R* changes sides one more time (fifth and 
ninth cross-flow planes) and thus a complete wave length 
of flow asymmetry is formed between the third and ninth 
cross-flow planes. Strong spatial shed vortex exists in the 
flowfield. This solution is strongly similar to the unsteady 
asymmetry local-conical flow solution at different time 
steps which is depicted in Fig. 8 on a cylinder with the 
axis of the cylinder representing time. The behavior of the 
flow asymmetry over one period in Fig. 8 is qualitatively 
similar to the behavior of the flow asymmetry over one 
wave length in Fig. 7. Figure 9 shows the total-pressure- 
loss contours and surface-pressure coefficient at different 
axial stations for the case of Fig. 7. The solutions at axial 
stations of X/L = 0.2 and 0.9 are almost the same (the 
total pressure losses are drawn to a scale given by the 
ratio of the circular diameters at X/L = 1 station and the 
local axial station). The flow asymmetry between these 
two stations represents a full wave length. 

Circular Cone-Cylinder Configurations 

To address the issue of the effect of cylindrical af- 
terbody length on the flow asymmetry a cylindrical af- 
terbody of different lengths is added to the unit-length 
conical forebody. The flow around the resulting cone- 
cylinder configurations is solved with the flow conditions 
of a = 40, Moc = 1.4 and R. = 4x 10 6 , which are the 
same flow conditions of the isolated unit-length cone of 
Fig. 2. The lengths of the cylindrical afterbody are cho- 
sen as 1, l.S and 2 and the corresponding grid sizes are 
taken as 161x81x65,161x81x69 and 161x81x72; re- 
spectively. The source of flow disturbance is the same 
short duration 2°-side-slip disturbance. The computed 
total-pressure loss for these cases are given in Figs. 10, 
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11 and 12; respectively. For the cone-cylinder configura- 
tion of 1:1 (cone length: cylinder length), Fig. 10 shows 
a very strong asymmetric flow on the cone, in compar- 
ison with the flow asymmetry of the isolated cone of 
Fig. 2, and on the cylindrical afterbody as well. It should 
be noted that inside the conical s h ock surrounding the 
cone-cylinder configuration, subsonic flow regions exist 
and hence the downstream cylindrical-afterbody bound- 
ary has an upstream effect. The cylindrical afterbody has 
dual effects which increases the flow asymmetry, the first 
is due to the cone -cylinder juncture and the second is due 
to the increase of the local angle of attack of the leeward 
side of the cylinder. Both of these effects increase the 
spatial growth of the flow asymmetry. 

For the cone-cylinder configuration of 1:1.5, Fig. 11 
shows stronger forebody asymmetry (first five cross-flow 
planes) in comparison with that of Fig. 10. On the cylin- 
drical afterbody, three cross-flow planes are only shown 
on a portion of its length. It should be noted here that 
slight flow unsteadiness has been detected during the 
computations. In Fig. 12, two snap shots are shown 
for the unsteady asymmetric flow solution of the cone- 
cylinder configuration of 1:2. Since the time step for a 
stable time-accurate solution of this case does not exceed 
10' 5 , it is computationally prohibitive to fully solve this 
case. 

Next, we show a comparison of the computed results 
with available experimental data. For this purpose, we 
consider the cone-cylinder configuration of 0.5:0.5 which 
was experimentally tested by Landrum*. The configura- 
tion angle of attack is 46.1°, the Mach number is 1.6 
and the Reynolds number based on the total configu- 
ration length (cone + cylinder) is 6.6 x10 s . The cone 
semi-apex angle is 9.5°. The problem is solved using a 
grid size of 161x81x65. Figure 13 shows the surface- 
pressure coefficient along with the experimental data, the 
total-pressure-loss contours and the total Mach-number 
contours at the axial stations of 0.075, 0.125, 0.225, 

0.475, 0.575 and 0.775. First, the computed and mea- 
sured surface-pressure coefficient are in good agreement 
on all the axial stations. Second, by studying the total- 
pressure-loss contours along with the total Mach-number 
contours, it is seen that flow asymmetry starts slightly at 
X/L = 0.0750, and spatially grows in the downstream di- 
rection. Moreover, the asymmetry changes sides in the 
downstream direction. This comparison conclusively val- 
idates our computed results and the grid size. 

Concluding Remarks 

The unsteady, compressible, thin-layer Navier-Stokes 
equations are used to obtain three-dimensional, asymmet- 
ric, vortex-flow solutions around cones and cone-cylinder 
configurations. Several important issues are addressed 
in the present study. By increasing the flow Reynolds 
number for flows around a cone, we have shown that 
the flow asymmetry becomes strong and changes sides in 


the downstream direction. For the high-ReynoIds flows, 
the spatial asymmetric flow develops in a wavy manner, 
which is qualitatively similar to the temporal asymmetric 
flow development of the locally-conical solutions, where 
the flow asymmetry develops in a periodic manner. By 
adding a cylindrical afterbody to the conical forebody, 
the flow asymmetry becomes stronger in comparison with 
that of the isolated cone. As the length of the cylindri- 
cal afterbody is increased, the flow asymmetry becomes 
stronger and unsteady. All these flow solutions have been 
obtained by using a short-duration side-slip disturbance. 
Finally, the computed results and grid used are conclu- 
sively validated. 
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Fieure 2. Asvmmetric flow solution around a cone of unit length, short-duration side slip. 
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Figure 8. 


Unsteady asymmetric locally-conical flow solution at different time steps within one f J 
cycle (cylinder axis is a time axis), At = 10" 3 . 4 ^ 
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Figure 10. Asymmetric flow solution around a cone-cylinder configuration 1:1. 
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-igure 11. Asymmetric flow solution around a cone-cylinder configurator 1:1.5. 
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Figure 12. 


Unsteady asymmetric flow solution around a cone-cylinder configuration 1:2, 
9767, At = 10" 5 . 
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